function model=update_alpha_1(model, n, MaxFun, data)


%options = optimset('LargeScale','on', 'MaxFunEvals', MaxFun, 'Display', 'off');

options = optimset('LargeScale','on','GradObj','on','Hessian','on', 'MaxFunEvals', MaxFun, 'Display', 'off');

x0      = model.alpha_1(n,:);
temp    = fminunc(@L_alpha_1, x0, options, model, n, data);

model.alpha_1(n,:) = temp;

end
